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ABSTRACT 

This  paper  is  an  attempt  to  provide  an  insight  into  the  rather 
complicated  problem  of  heading  error  in  spin  resonance  magnetometers. 
The  particular  system  investigated  is  the  self  oscillating  rubidium  vapor 
magnetometer,     Theory  is  developed  which  shows  the  magnitude  of 
heading  error  to  be  a  function  of  many  parameters;  however,  only  two 
of  these  parameters,  the  optical  pumping  rate  and  the  r-f  feedback, 
are  singled  out  for  experimental  varification.    The  results  are  encouraging 
in  that  they  show  that  the  magnitude  of  heading  error  follows  a  definite 
pattern;  consequently,  it  should  be  possible  to  introduce  compensating 
networks  to  reduce  or  possibly  eliminate  entirely  the  orientation  effects. 

The  material  for  this  paper  was  gathered  while  the  author  was  a 
guest  of  Varian  Associates  in  Palo  Alto,  California.    Thanks  are  expressed 
to  the  staff  of  the  instrument  division,  and  especially  to  Dr»  J.  A,  Arnold; 
for  their  suggestions  and  technical  help. 

The  writer  also  wishes  to  thank  Professor  Carl  E„  Menneken  of  the 
U.S.  Naval  Postgraduate  School  for  his  constructive  criticisms  in  the 
preparation  of  this  paper. 
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I.    INTRODUCTION 

Precision  magnetic  field  measurements  are  being  used  in  an 
increasing  number  of  military,  scientific,  and  industrial  applications. 
Many  of  these  applications  are  such  that  the  measurements  must  be 
obtained  from  a  moving  platform.    This  means  that  the  instrument  used 
must  be  compact;  it  must  be  rugged  enough  to  sustain  forces  such  as 
those  encountered  when  launching  or  recovering  satellites;  it  must  use 
little  power;  and  it  must  be  more  sensitive  than  existing  operational 
devices . 

The  recently  developed  rubidium  vapor,  spin-resonance,  magne- 
tometer meets  these  requirements.    It  has  already  been  flown  as  part 
of  the  instrument  package  in  several  satellites  and  its  performance  met 
or  exceeded  expectations.    It  is  anticipated  that  the  magnetometer  will 
perform  equally  well  in  the  military  applications  for  which  it  is  presently 
being  packaged . 

Like  most  new  devices,  there  are  engineering  problems  that  must 
be  solved  before  the  magnetometer  is  completely  operational.    The  one 
with  which  this  paper  is  concerned  is  heading  error-  the  variation  of 
field  strength  reading  as  the  orientation  of  the  magnetometer  is  changed. 
It  will  be  shown  that  the  orientation  effects  are  fundamental  to  spin- 
resonance  magnetometers  using  an  absorption  media  having  nuclear  spin 
other  than  zero.    It  will  also  be  shown  that,  to  some  extent,  the  magni- 
tude of  the  heading  error  is  a  function  of  design  parameters  which  may  be 
altered  to  suit  a  specific  application  of  the  instrument. 


II.    THEORY 
The  particular  magnetometer  to  be  discussed  is  an  oscillator 
whose  frequency  is  determined  by  the  energy  difference  between  the 
Zeeman  split,  hyperfine  levels  of  the  5s  state  of  the  rubidium  85  atom. 
This  energy  difference  is  in  turn  determined  by  the  magnitude  of  an 
external  magnetic  field,  in  this  case  the  earth's  field.    By  utilizing 
appropriate  techniques  to  measure  the  oscillator  frequency,  absolute 
field  measurements  to  a  few  gamma,  or  field  perturbations  to  ,01  gamma 
may  be  made.    To  describe  the  operation,  it  will  be  necessary  to  investi- 
gate the  techniques  used  to  perform  the  basic  operating  steps  of  aligning 
and  tipping  the  nuclear  magnetic  moments ,  and  then  detecting  and 
measuring  the  resulting  precession  frequency. 

A,    Alignment  of  the  Nuclear  Magnetic  Moments 

The  alignment  of  the  magnetic  moments  is  accomplished  by  utilizing 

I 
the  optical  pumping  techniques  described  by  Kastler    .    A  pyrex  cell 

containing  rubidium  vapor  is  irradiated  with  the  light  from  an  excited 

rubidium  lamp      Normally  this  stimulates  transitions  out  of  the  S  to  the 

P  states  corresponding  to  the  quantum  mechanical  restriction    /\m  =  ±  I 

By  using  circularily  polarized  light,  however,  the  transitions  can  be 

restricted  to  either  4mz  +  l    or/\m=-|      depending  on  the  sense  of  circular 

2 

polarization    .    For  light  polarized  such  that  Am=+i,  the  transition 

A.  Kastler,  Proc.  Phys .  Soc .  ,  A67,  863  (1954) 

2 
G .  Herzberg ,  Atomic  Spectra  and  Atomic  Structure 


probabilities  out  of  the  various  sublevels  of  the  ground  state  have  been 
calculated  and  are  tabulated  in  Fig.   1.    Atoms  in  excited  states  return 
almost  immediately  to  the  ground  state  by  emission  of  a  photon.    Since 
the  return  transition  is  not  restricted  toAm:+l  «   and  since  there  is 
considerable  mixing  in  the  excited  sPy    and    P3,  levels,  (a  condition 
that  may  be  enhanced  by  the  addition  of  a  buffer  gas) ,  the  return  to 
the  ground  state  is  random.    It  is  apparent  then  that  at  some  time,  T    , 
the  upper  levels  of  F=3  and,  to  a  smaller  extent,  the  lower  levels  of 
F=2  will  be  over  populated,  while  the  lower  levels  of  F=3  and  upper 
levels  of  F=2  will  be  nearly  depleted. 
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FIG.  I 
TRANSITION    PROBABILITIES  OUT  OF  TH£  £Sy&    TO    THE    zP</e    STATE 
hi ITH  THE  QUANTUM    MECHANICAL    RESTRICTION     Am  =  +  I 


The  optical  pumping  process  is  illustrated  in    Fig's  2  and  3.    A 
model  having  zero  nuclear  spin  is  used  for  the  illustrations  in  order  to 
avoid  the  confusing  complexity  of  the  actual  R,    85  energy  level  diagram, 
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however,  the  principles  involved  are  basically  the  same.    Fig.  4  is 
the  actual  energy  level  diagram  of  the  ground  state  and  shows  the  field 
dependent  radio  frequency  transitions. 

The  time  necessary  to  complete  the  pumping  process,  Tp ,  will  be 
dependent  not  only  on  the  intensity  of  the  light  source,  but  also  on  the 
degree  to  which  the  light,  as  viewed  parallel  to  the  earth's  field      H0    > 
is  elliptically  polarized.    Right  hand  circularily  polarized  light  may  be 
described  by  the  equations: 

(D        Px=  A  cos  oji 
Py  -  A  sin  oot 
Referring  to  Fig,  5,  if  the  circularily  polarized  source  is  rotated  with 
respect  to    H0     ,  the  observer  now  sees  elliptically  polarized  light. 

^e/USMr 

/  fig.  5 

*  rotation  of  light  axis 

with  respect  to  h0 

■jr 

-Z 

Thus,  for  any  angle,   •&-    c  the  polarization  as  seen  by  the  observer  is 
described  by 

(2)     PA  z  A  cos-e-cos  coi 
Py  =  A  sin  u>t 

From  (2)  it  is  apparent  that  the  degree  of  light  polarization  has  a  cos-e- 

variation,  and  for   -0-=9O°,  the  polarization  is  linear  and  no  optical 

pumping  is  achieved.    Since  the  signal  to  noise  ratio  of  the  output  is 

dependent  on  the  extent  to  which  the  sample  is  pumped,  the  sensitivity 


will  also  have  a  cos  -&■   dependence.    The  plane,   -6- =90    ,  is  referred 
to  as  the  equatorial  dead  zone. 

Note  that  had    He     been  rotated  and  the  light  orientation  held 
fixed,  the  results  would  have  been  identical. 


B.    Spin  Displacement 

The  result  thus  far  has  been  to  produce  a  net  alignment  of  the 
individual  nuclear  magnetic  moments,  /cF.    ,  at  an  angle  to     Mo 
dictated  by  quantum  energy  considerations.    The  magnetic  vectors 
associated  with^r.  precess  around    H0    xn  an  incoherent  manner 
such  that  if  a  total  magnetization  vector,     M     ,  is  defined  where 


VM  =  ^  ^  f£        J  then 


(3)     M  =  Mx  J   +  My  J  +  Mz  k 

and  in  this  case    Mx~  My  =  0       .    The  purpose  of  this  section  is  to 
show  how    M     is  tipped  with  respect  to   H0  ,  and  subsequently 
precesses  about   r\0   at  the  Larmor  precession  frequency, 


■LOCUS  OF  /if i 


(a)  M    J        (b) 

FIG.  6 
(do)  LOCUS  OFAfI  before  optical  pumping  (neglecting  the  small  effect 

OF    30LTZMAN    DISTRIBUTION)        (b)  ALIGNMENT    OF   7tfC    AND   RESULTANT 
ft     DUE   TO    OPTICAL   PUMPING 

If    M      is  initially  assumed  to  be  displaced  from    H0   by  an 


dt 
with  the  resulting  rotation  in  a  plane  perpendicular  to  both  M  and    H0 


angle,     c*     ,  a  torque  given  by  - — =«fMxH0     will  be  produced, 

dt 


_      3 

Now  assume  a  coordinate  system  rotating  with  angular  velocity  -co    . 

With  respect  to  this  new  coordinate  system 

(4)         iB=^M+(SxM) 
d  "t  at 

»("lMxH.)  +  (  i5  x  Wm) 


where  f-/.~  <±L      is  the  effective  field  and  is  stationary  when  viewed  in  the 
i 

rotating  frame . 

Thus  the  requirement  that    M     be  displaced  from    Ho  is  clearly 
illustrated,  and  a  method  of  maintaining  this  displacement  against 
the  pumping  action  of  the  rubidium  lamp  is  suggested.    Let  a  second 
field,    H,     ,  which  rotates  with  angular  velocity  -  6U  ,  be  superimposed 
on  the  system.    As  viewed  in  the  rotating  frame,  hi,      will  be  stationary 
and  the  effective  field  becomes 

(5)      Heff  =(/70-  £)+   H, 

If  the  coordinates  are  chosen  such  that 


(6)    H,ffm    (Ho-^J  k  +  h,  7 


then  SM     vanishes  only  for  o<  =90    ,  consequently    s!M     cannot  vanish, 

it  dt 

and   ^    must  remain  displaced  with  respect  to     Ho   •    At  resonance, 

600  =  ^Ho     f  Hocc  s  Hit     /  and-T0-  —   which  is  the  Larmor  precession 
frequency . 


^I.I.  Rabi;  W.  F.  Ramsey;  and  J.  Schwinger,  Reviews  of  Modern 
Physics,  167,  Vol  26  (1954) 
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The  tipping  of  M    represents  a  change  in  energy  ,^H   /  and  since 

this  is  restricted  to  a  quantum  change,,  it  remains  to  be  shown  that    f0 

is  the  frequency  corresponding  to  the  transitions  discussed  in  section  A. 

Neglecting  second  order  and  higher  effects ,  and  referring  to  the  energy 

level  diagram  of  the  2Sy     state,  (Fig.  4): 

(?)    w(f  k)  -  N  -  H0  A  Where /c^      is  the  component  of  /*F 

^H  in  the  field  direction 


(8) 


Kl 


(F^i)-  N(rms)m-H0    -i^(mi.^j) 


and 

since 

\mi 

-^;l 

=  1 

(9) 

1*1  r 

2TT 

- 

.  4GG  NQ  (iou)      cycles/ sec 

In  actuality,  the  requirement  that    H(    be  a  field  rotating  at 
angular  velocity  -co  is  satisfied  by  a  plane  polarized  wave,  since 
a  plane  polarized  wave  can  be  considered  to  be  the  combination  of  two 
circularily  polarized  waves  rotating  in  opposite  directions.    Also,  H, 
has  thus  far  been  assumed  to  be  in  a  plane  perpendicular  to     Ho   /  or 
in  the  notation  of  the  previous  section,   -G-    =90    .     For     -0-^90° ,  only 
the  component  of    h\\    which  is  orthogonal  to    H0    is  useful  in  inducing 


transitions.    The  useful  field,  then,  is  given  by  Mj=   Ht  s(n&  ,  and  the 
overall  sensitivity  of  the  magnetometer  will  have  a  sin  ■©■     as  well  as 
acos^      orientation  dependence. 


10 


C.    Detection  of  the  Precession  Frequency 


If  the  concept  of  a  rotating  set  of  coordinates  is  dropped,  the 
rubidium  vapor  sample  can  be  characterized  as  a  total  magnetization 
vector,    ffi    ,  precessing  around    Ha     at  the  Larmor  frequency.    The 
problem  now  arises  as  to  how  to  detect  this  frequency.    The  most 
successful  method  thus  far  has  been  a  cross  beam  apparatus  func- 
tionally equivalent  to  that  shown  in  Fig .  8 . 
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In  a  manner  analogous  to  that  of  section  A,  it  can  be  shown  that 
a  second  circularly  polarized  light  source  with  its  axis  in  the  plane  of 
H ,    will  also  have  an  effect  on  the  rubidium  sample.    The  sample  will, 
in  fact,  be  optically  pumped  by  the  second  lamp;  however,  the  transition 
probabilities  describing  the  pumping  action  will  have  oscillating  co- 
efficients, thus  the  light  reaching  the  photo  detector  will  be  modulated 

at  the  precession  frequency.    For  the  two  energy  level  model,  Bell  and 

3 
Bloom     have  derived  the  following  phenomenological  equations  for  the 

detector  output: 

3W.E.  Bell  and  A.  L.  Bloom;  Phy.  Rev.  ,   107,   1560  (1957) 

11 


,2-,/p 


<10>     5X=   VHMr5a  OfeAa>rj  z     cos(cot  +  >,>) 

H-(SzAoo)£+iaH?S,Sa 


(H)      S£  =      Mrlu-  fS2Aco)  ] 


*-,% 


/+  (SaAGjf+1*H*S,S£ 
where 

'i -(*"*♦£'" 

P  2  is  the  rate  of  optical  pumping  parallel  to    H0 

P^  is  the  rate  of  optical  pumping  perpendicular  to   /-/< 

77  is  the  thermal  relaxation  time 

T2  is  the  transverse  relaxation  time 

My  is  the  equilibrium  population  difference 

At  resonance  the  signal  becomes 


(12)       ^  a      jHlS2Mr        CC)S^t  +  yj 


The  system  actually  used  in  the  magnetometer  does  not  use  a 
separate  monitoring  lamp  but  uses  a  single  lamp  offset  from   Ho  so  that 
the  beam  has  components  in  both  the  X  and  Z  directions.    This  means 
that  the  monitoring  light,  and  consequently  the  signal  intensity,  will 
have  an  additional  sin-©-     orientation  dependence  that  must  be  combined 
with  the  factors  previously  discussed.    The  overall  magnetometer 
sensitivity,  as  related  to  orientation,  is  given  by  • 


12 


Fig.  9  is  a  plot  of   (13)  for  K=1.5. 

To  make  the  system  a  self  oscillator  the  output  of  the  photo 
detector  is  amplified  and  fed  back  to  the  r-f  coil  to  produce  the  H, 
field . 

The  complete  magnetometer  is  shown  in  Fig.  10. 
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CROSS-BEAM    RUBIDIUM  VAPOR     MAGNETOMETER 
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Ill .     HEADING  ERROR 
A.    Definition  of  Heading  Error 

In  this  thesis  the  term,  heading  error,  is  defined  as  the  difference 
between  the  oscillator  frequency  when  the  magnetometer  is  oriented  at  an 
angle,  <p  -  0  ,  and  when  it  is  oriented  at  a  second  angle,  <p  -  <$  .  The 
zero  reference  is  the  polar  dead  zone,  <±  -  2/o°  •  Thus  $>-4>  -  ■&■  ,  the 
offset  from  Ho  >  however,  -Q-  now  carries  a  sign  which  eliminates 
ambiguities  as  to  the  direction  of  offset.  The  axis  of  rotation  is  always 
the  East-West,  or  latitude,  line. 
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B.     Maximum  Heading  Error 

At  present  there  is  no  analytic  expression  from  which  the  exact 
value  of  heading  error  may  be  determined.  Fo  a  given  instrument  the 
limits,  or  maximum  heading  error,  however,  may  be  calculated  from  the 

frequency  spectrum  obtained  by  a  straight  forward  application  of  the 

4 
Breit-Rabi  equation 


14)  wM(F  ,  =  -JUL--  4i  Hm  1  4*    /,4 


'(81+ 1) 


4m    x+><2 


£T+l 


where     AN  -  2TT  f\  Af  and      X  =    (  "Aj/jr  +  Ar/r  )  H 

AlN 

An  approximate  plot  of  (14)  is  shown  in  Fig.   12. 


For  small  values  of  x,  to  second  order 
(15) 


FIG, IS 

VARIATION  OF  ENERGY 
NITH  MAGNETIC  FIELD 


f: 


m-X+X£    z     I  -r   JJ2-(-Aj/T  4  Aj/r)— , 
2T+\  21  + 1    I  J  JJ  AN 


<£-' 


(21+1) 


•K-**-*f& 


G.  Brelt  and  I.  Rabi ,  Phys  .  Rev.  ,  38,   208  2  (19  31) 
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Substituting  (15)  back  into  (14) 


(16)  ifikfafiA*  --AbL AjHm  ±  m[i+  ?m   t-Aj/ +  Aj/\\±. 

[a    (cT+,)zJi     -*-x]Xn2  J 


/4  k/ 
and 


(17)    Kw  fam,)  -  NMft  ntyn-df.  H(mr^){/^  }  y  ^LMirma) 

Using  (£r+,;2  ^W 

Then 


(18)  J^  (»i,-hia)Ho  =  '^^Ho  *  W2> 

(19)  f0  -    .  4££   ^f/o6)   *      ^^   kc/s 

Note  that  (19)  is  identical  to  the  first  order  approximation  of  Part  II. 
Looking  at  the  second  order  effects: 

(21)  (-A&  +/'fif   «"  (m«  _  „£)  =  ,f,ah  (m*-m?)=  h^ 

where 

^m"mn   =    ^3  ,J'lr^j"S     for    F*3 

(22)  r    3,lrlrJ      "for    F^£ 

combining  (20),   (21),  and  (22)  the  resulting  frequency  spectrum  is  that 
shown  in  Fig .   13. 
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FIG. IS 
RESONANT   FREQUENCY    SPECTRUM    OF  Rb  8S 

Fig.   13  shows  that  the  R^  85  resonant  frequency  spectrum  is 
composed  of  not  one  but  ten  distinct  frequencies.    The  four  lines  cor- 
responding to  the  F=2  transitions  are  of  little  consequence,  however, 
as  it  can  be  shown  that  the  intensity  of  these  lines  is  small  compared 
to  those  associated  with  the  F=3  transitions. 

The  frequency  of  oscillation,  then,  can  be  expected  to  be  some 
ensemble  average  of  the  F=3  spectrum.    If  Fig.   13  is  now  redrawn  to 
show  the  relative  line  intensities  when  the  rubidium  vapor  is  polarized 
withAmr+i     andAm =- 1  transitions ,  (Fig.   14),  it  is  apparent  that  the 
frequency  of  oscillation  will  be  different  depending  on  the  sense  of 
rotation  of  the  circularily  polarized  light. 


kAf  > 


Af  4 


Am=+I 


Am  =  -I 

FIG.  IR- 
RELATIVE   INTENSITIES   OF  THE  SPECTRUM 
LINES 


A^n=  +  I    (broa.denad  line) 
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LetAf    be  the  difference  in  the  arithmetic  mean  frequency      f0     , 

and  the  frequency  of  the  resonance  line  maximum.     From  Fig.   14  it  can 

be  seen  that  under  conditions  of  complete  optical  pumping  the  frequency 

spectrum  for  F=3  would  consist  of  only  one  line,  f=f  -89.    Since  Af 

o 

is,  in  effect,  a  measure  of  the  competition  between  the  optical  pumping 
and  the  magnitude  of  the  effective  r-f  feedback,  it  corresponds  to  the 
previous  definition  of  heading  error.    The  maximum  heading  error  obtain- 
able, therefore,  is  approximately  90  cycles/second.    It  should  be  noted 
that  this  value  is  considerably  greater  than  that  physically  obtainable 
with  this  system. 
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C .    Theoretical  Values  of  Heading  Error  Versus  Orientation 

The  consequences  of  Af    would  not  be  serious  if  it  had  one  value, 

Af     ,  throughout  one  hemisphere,  and  changed,  as  the  magnetometer  is 
rotated  through  the  equatorial  dead  zone,  to  a  second  constant  value, 

A-J^    .    This  type  of  error  could  easily  be  corrected  by  a  minor  adjustment 
of  the  readout  equipment.    Unfortunately,  due  to  line  broadening  effects, 
Af   is  not  constant  for  any  appreciable  range  of  -Q-  .    Some  of  these  effects 
have  been  treated  analytically;  however,  at  the  present  time,  these 
treatments  are  for  single  lines  and  do  not  apply  to  the  rotating  system 
under  consideration. 

It  will  suffice,  as  far  as  the  scope  of  this  paper  is  concerned,  to 
consider  all  of  the  line  broadening  effects,  except  the  optical  pumping 
rate  and  the  r-f  feedback,  as  part  of  the  transverse  and  longitudinal 
relaxation  time  constants,  T     and  T„;  and  then  look  at  the  resonant 
line  shape  as  a  function  of  H,  /  Px  ,  Pz  ,  and  ■&■  .    It  should  be  kept  in 
mind  that  anything  that  broadens  the  resonance  line  also  brings  the 
oscillating  point  nearer  to  f    and  consequently  reduces  heading  error. 
Assuming  H     to  be  a  rotating  field,  the  complete  solution  of  the 

Bloch  equations  for  the  transverse  component  of  the  magnetization 

5 
vector,  when  corrected  to  include  optical  pumping,  is 

(23)       Mx=-i-MTV«52     Sa(<*>ru>)2M,cos<*)t    +  3H,smooi 
2  i+S*(u0-cj)£-h  fH*S,Ss 


J.  A.  Pople ,  W.  G.  Schneider,  and  H.J.  Bernstein,  High-Resolution 
Nuclear  Magnetic  Resonance 
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where      S|  =  />„  +  Pz  +  J_\~'    a,nc|      S£  =    (P^  +  Pt+±  \~' 

The  effect  of  H    being  a  linear  rather  than  an  oscillating  field  does 

not  effect  the  validity  of  the  solution;  however,  the  solution  now  contains 

both  in-phase  and  out-of-phase  components  with  respect  to  the  applied 

field,  H,=  sHTC0S-&-  •    These  components  can  be  described  in  terms  of 

the  Bloch  susceptibilities ,  %     and  %       ,  where  %    is  the  in-phase 

component. 

<24>   My,*  &x'H,cosoA, 

and   Aj      is  the  out-of-phase  component 

(25)  Mx2=2%"Hl<sihu>h 

Since   MT  -  %0  K,        and    jfc.-  X'+J  X* 

(26)  #'  =    ±  Mr       $?(coo-00) 

2 


l+%'(cu0-cof+1iHt*$lS. 


(27)  v"=i    m        § 


a 


The  mean  rate  of  absorption  per  unit  volume  is  given  by 


(28)       A     =_M      dHx/dt 


3 


/  +  S*  (co.-co)2+  i*u,ss,sc 
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Since  the    coscot  bin  <-*>i       term  averages  to  zero,  and  the  sin    cot 
term  averages  to  1/2.    This  shows  that  power  is  absorbed  only  by  the 
out-of-phase  ,  )£     ,  component.    Since  the  magnetometer  operates  on 
the  absorption  mode,  the  signal  intensity  is  proportional  to 


<29)       Mxp   =       MT  iH,Sa 


'+SJ7a,fl-w)a+  l'H,*S,Sa 


sm  cot 


This  can  be  obtained  from  equation  (10)  given  in  section  II;  how- 
ever, now  it  is  clear  that  a  90  degree  phase  shift  is  necessary  in  a  self 
oscillating  system,  and  it  is  possible  to  introduce  the  dispersion  mode 
which  comes  from  the  in-phase,  or    %     component. 


d00 

FIG. 15 
ABSORPTION    LINE 


Fig.   15  is  a  plot  of  the  absorption  curve.    The  Lorentzian  shape 

does  not  hold  for  the  rubidium  resonant  line  since  equation  (30)  is  valid 

for  transitions  between  two  lines  only,  however,  anything  which  broadens 

this  line  will  also  broaden  the  rubidium  line. 

Using  the  approximation,    T     =  T       equation  (29)  reduces  to 

1         ^ 
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(30)    My_  =        Mr-/^.  S 


xg  = 


It  is  easily  shown  from  (30)  that  the  maximum  value  of    Mxe     occurs  for 

H,=   t4-  ,  and  that   My  =  0    for  ^,=  0     or    tf,  =  *0    .    The  lind  width  is 

proportional  to  H    .    Since     H.  =  Hr  Sih-Q-         ,  and  since  it  has  been 
1  ' 

postulated  that  heading  error  is  inversely  proportional  to  line  width; 

it  can  be  further  postulated  that  heading  error  is  a  function  of 
{'-lAsmolj    for    W,>^     ■ 

The  effects  of  S(£)  variations  should  not  be  as  severe  as  those  of 
H(G)  since  S=(£+/^+^  )        =     f /£(»>*>)  COS e  +  B(mM)*in-e-  +  -=j=  1 
and  for     -i-    /<    R.  x        ,  the  range  of  S  would  beiL=  l  ,{t  J  for  &an°4saCinQ 
The  width  of  the  resonance  line  is  inversely  proportional  to   v$      ; 
therefore,  if  no  AGC  is  used  the  heading  error  should  also  be  a  function 
of    //  -J6sinBe\J  . 

The  heading  error  versus  orientation  curve  for  the  absorption 
oscillation,  then,  can  be  expected  to  be  similar  to  Fig.   16.     Fig.    17 

is  a  plot  of  the  dispersion  curve.    By  utilizing  the  appropriate  phase 

o 
shift,  0  or  180    ,  the  magnetometer  can  be  made  to  oscillate  on  either 

of  the  two  peaks.    These  peaks  now  occur  at  a  frequency,    go     ,  for  which 

H  -  yi+  S2  (ui-cJo)2"      -v   1 OJ  -  cJo /  .    Thus  the  frequency  separation  of  the 

'  IS  1 

maxima  is  directly  proportional  to  the  magnitude  of  H    .    As  H     is 

1  1 

decreased  the  dispersion  oscillation  frequency  does  not  reduce  to  f 

o 

but  to  the  absorption  frequency.    This  means  that  the  heading  error    will 
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POSSIBLE     HEADING    ERROR    VERSUS    ORIENTATION    FOR 
THS    ABSORPTION      OSCILLATION 
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not  be  symmetrical  about  the  equatorial  dead  zone.    There  may  even 
be  an  optimum  choice  of  /-/rsuch  that  as  the  instrument  orientation  is 
varied  throughout  one  hemisphere,  the  effects  of  W(  on  line  broadening 
and  the  dispersion  peaks  separation  cancel  and  the  heading  error  will 
be  essentially  constant.    A  heading  error  versus  orientation  plot  could 
conceivably  be  that  shown  in  Fig.  18. 
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POSSIBLE   HEADING  ETRROR    VERSUS   ORIENTATION 
FOR    THE    DISPERSION     OSCILLATION 
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IV.     EXPERIMENTS 

A.    Polarizer  Rotation 

o 
In  the  experiments  to  follow,  a  180     rotation  of  the  magnetometer 

was  simulated  by  changing  the  polarity  of  the  pumping  light.    This  has 

several  advantages:     (1)  The  effects  of  magnetic  components  in  the 

system  are  minimized.     (2)  The  light  polarity  can  be  changed  in  a 

fraction  of  the  time  necessary  to  rotate  the  entire  system,  thus  the 

effect  of  field  variations  is  minimized.     (3)  The  system  components  can 

be  assembled  in  more  accessible  locations.     (4)  The  necessity  to  manually 

switch  the  r-f  connections  when  going  from  one  hemisphere  to  another  is 

eliminated . 

The  circular  polarizer  is  composed  of  a  linear  polarizer  followed  by 

a  quarter-  wave  plate.    The  polarizer  and  quarter-wave  plate  are  aligned 

such  the  quarter-wave  plate  sees  two  equal  components  of  linearly  polarized 

o 
light,  one  of  which  is  passed  unaltered,  the  other  delayed  90    .    The 

o 
polarity  of  circular  polarization  is  changed  by  a  90    rotation  of  either  the 

linear  polarizer  or  the  quarter-wave  plate. 

The  reason  that  the  sense  of  rotation  of   f-f,    /  does  not  have  to  be 

reversed  in  order  to  maintain  an  oscillation  when  the  light  polarity  is 

reversed  is  easily  explained  with  reference  to  the  classical  vector  model. 

Fig.   19a  shows  the  instantaneous  position  of  the  angular  momentum 

vector,    F     ,  when  2  W,  cos  cot  =  2^    .    For  4m=+|  ,    ~pp      is  oppositely  aligned 

with  respect  to  F     ;  thus  the  precession  and  tipping  torques  given  by 

£,  =  >F  *  K>  and    L£  -/iF  X  H,       ,  result. 
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(b) 

riG.ia 

PHASE    RELATIONSHIPS    IN  THE    SELF- OSCILLATING   MAGNETOMETER 


Fig.   19b  shows    f/0    reversed.    Now  in  order  to  produce  a  tipping 
torque  and  still  maintain  the  9  0     phase  lag  necessary  for  oscillation, 
it  is  evident  that  the  polarity  of   H,      must  be  reversed. 

In  Fig.  19c,  the  polarity  of  light  is  reversed  instead  of    HQ     . 
Optical  pumping  now  is  in  the  negative  energy  direction  and  tends  to 
align  /tF      with   H0    ;  however,    /rF     is  also  aligned  with    F      since 

Atn~-I   i  and  the  net  result  is  a  reversal  of  precession  direction.    The 
conditions  for  oscillation,  then,  are  maintained  without  changing  the 
polarity  of     Ht     . 
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B.    Experimental  System 

The  components  illustrated  in  Fig.  10  were  enclosed  in  a  cylindrical, 
aluminum  container,  a  section  of  which  was  removed  to  allow  free  access 
to  the  components  and  to  facilitate  rotation  of  the  quarter-wave  plate . 

The  temperature  of  the  gas  cell  was  maintained  at  approximately 
40    C  by  a  60  cycle  heater  which  consisted  of  a  bifilar  winding  of  5  turns 
of  nichrome  wire  wound  over  the  r-f  coil.    A  calibrated  thermistor  was 
used  to  moniter  the  temperature. 

The  complete  system  mounted  on  a  wooden  rotatable  platform  is 
shown  in  Fig.  20. 


FIG.  20 
EXPERIMENTAL    APPARATUS 
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C .    Line  Shape 


4W££P    CONTROL 


fig.  ai 

EXPERI MENTAL   <S?T  UP    FOR  OBTAINING,    LINE    SHAP? 


The  driven  system  diagramed  in  Fig.  21  was  used  to  obtain  the 
resonant  line  shape.    Unlike  a  self-oscillator,  the  output  of  a  driven 
system  is  neither  the  dispersion  nor  the  absorption  line,  but  a  combination 
of  the  two.    This  is  possible  since  the  phase  angle  between  the  externally 
applied  r-f  field  and  the  detector  output  is  not  fixed,  but  varies  with  fre- 
quency such  that  the  detector  output  is  a  maximum.    Although  this  pre- 
vents desirable  absorption  band  width  measurements,  it  does  give  useful 
information  in  that  the  absorption  line  is  predominant  when  the  r-f  level 
is  below  the  saturation  value  while  the  dispersion  line  is  predominant 
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for  values  of  r-f  above  the  saturation  level. 

A  15  to  20  second  linear  sweep  was  used  to  drive  both  the  scope 

horizontal  input  and  the  varicap  controlled  oscillator.    The  sweep  was 

generated  manually,  using  a  line  and  pully  arrangement  attached  to  a 

ten  turn  helipot,  R    .    The  second  pot,  R    ,  was  adjusted  to  obtain  the 

1  2 

desired  frequency  sweep. 

The  line  shapes  for  various  values  of  r-f  drive  were  photographed 
using  a  polaroid  camera  set  for  a  time  exposure.     Frequency  calibration 
of  the  photographs  was  obtained  by  disconnecting  the  scope  vertical 
input,  moving  the  remaining  spot  to  significant  positions  on  the  scope 
with  R    ,  and  making  multiple  exposures  to  record  these  positions  on  the 
photograph  while  recording  the  frequency  corresponding  to  each  spot 
position. 

Representative  results  are  shown  in  Fig.  22. 
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FIG.  2£ 
RESONANCE   LINE  SHAPE  AS  A   FUNCTION  OF 
R-F  DRIVE 
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D.    Heading  Error  Versus  Phase  Angle 

A  determined  effort  was  made  to  measure  the  phase  angle  between 
the  output  of  the  detector  and  the  input  to  the  r-f  coil.    Since  the  photo- 
cell gives  the  90  phase  lag  necessary  for  the  absortion  oscillation,  and 
since  the  ratio  of  the  coil  impedance  to  the  amplifier  output  impedence 
is  on  the  order  of  1/50;  this  measurement,  added  to  90,  would  have  been 
a  very  good  approximation  of  the  phase  angle  between  the  light  intensity 
and  the  r-f  field,  HI. 

A  comparator  capable  of  measuring  phase  shifts  of  approximately 
t  5°         was  designed  and  constructed.    Although  bench  tests  of  the 
device  were  satisfactory,  it  would  not  work  in  the  experimental  magne- 
tometer system.    In  the  light  of  later  developments,  it  appeared  that  the 
most  probable  reason  for  its  failure  is  that  an  external  phase  shifter, 
used  to  obtain  a  good  strong  oscillation,  produced  an  additional  90  shift 
which  was  completely  beyond  the  range  of  the  comparator.    The  signal 
thought  to  be  the  absortion  oscillation  was  subsequently  identified  as 
a  dispersion  oscillation 

The  system  finally  used  was  the  self-oscillator  shown  in  Tig.  23. 

The  desired  mode  of  oscillation  was  selected  with  the  0-180     active 

o 
phase  shifter,  and  the  phase  was  varied  approximately  13    either  side 

of  the  oscillating  point  with  a  series  L-C  network.    The  amplitude  of 

the  r-f  input  was  held  to  5mv  by  varying  the  output  level  of  the  active 

phase  shifter  to  compensate  for  the  drop  across  the  L-C  circuit. 
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RF    VOLTMETER 

^IS'VARICAP 
PHASE    SH/FTER 

FIG-.  £3 
EXPERIMENTAL    SETUP    FOP     OBTAINING    PHASE 
ANCLE    AND     ORIENTATION    EFFECTS 


FIG.  B.-4 
VARICAP  PHASE    SHIFTER   CALIBRATION  CURVE 
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FIG.  B5 
PHASE   SHIFT  VERSUS    BAf 
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E.    Heading  Error  Versus  Orientation 

Although  the  present  magnetometers  are  not  designed  to  oscillate 
at  the  dispersion  peak,  it  was  found  that  by  inserting  an  additional  90° 

0 

phase  to  either  cancel  the  90  introduced  by  the  photocell  detector  or 

e 
increase  it  to  180,  an  oscillation  was  obtained.    To  verify  that  this  was 

the  dispersion  mode     the  following  experiment  was  performed: 

(1)  With  the  magnetometer  initially  oscillating  at  a  suspected  dis- 
persion peak,  the  r-f  drive  was  reduced  until  the  oscillation  ceased. 
Depending  on  the  magnetometer  orientation,  the  minimum  r-f  level 
necessary  for  oscillation  was  5  to  10  mv.    The  r-f  drive  was  then  reduced 
still  further,  to  approximately  lmv. 

(2)  The  phase  was  then  varied  until  a  weak  oscillation  was  again 
obtained.    The  phase  shift  was  approximately-90   ,  thus  it  was  assumed 
that  this  was  an  absorption  oscillation.    To  check,  the  phase  was  held 
constant  and  the  r-f  level  gradually  increased.    The  signal  amplitude 
increased  to  a  saturation  point  and  then  decreased.    The  r-f  level  was 
further  increased  until  the  oscillation  either  completely  disappeared  or, 

in  some  cases,  reduced  to  a  violent  squegging. 

o 

(3)  The  r-f  was  held  constant  and  the  phase  angle  shifted  +90, 

to  its  original  value.    A  strong  oscillation  about  200  cycles/sec  higher 
than  the  absortion  oscillation  frequency  was  observed. 

Since  the  results  of  (1) ,   (2) ,  and  (3)  coincide  with  those  predicted 
by  the  Bloch  equations,  it  was  concluded  that  it  was  possible  to  obtain 
both  dispersion  and  absortion  mode  oscillations . 
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Data  was  obtained  to  determine  the  heading  error  for  both  modes  of 
oscillation.    The  experimental  system  and  procedure  was  the  same  as  that 
used  to  obtain  the  heading  error-phase  relationship,  except  that  the 
magnetometer  was  rotated  and  the  phase  held  constant.    The  frequency 
difference  observed  by  rotating  the  circular  polarizer  is  taken  as  £Afj  £e. 
it  is  assumed  that  heading  error  is  symmetric  about  the  equatorial  dead 
zone. 

The  results  of  this  experiment  are  tabulated  in  Fig.  26.    No  in- 
ferences should  be  drawn  from  direct  comparisons  of  oscillating  frequencies 
at  different  magnetometer  orientations  .    The  magnetic  field  gradient  in  the 
room  in  which  this  experiment  was  conducted  was  so  severe  that  it  was 
found  necessary  to  use  a  small  compensating  field  to  obtain  a  satisfactory 
absorption  oscillation,  and  the  amount  of  compensation  was  not  necessarily 
the  same  for  all  orientations.    Care  was  taken,  however,  to  insure  that  the 
total  field  vector  at  the  absorption  cell  remained  fixed  at 
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FIG.  26 
TABULATED    VALUES  OF  Af   VERSUS  ORIENTATION 
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V.     EVALUATION  OF  RESULTS 
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FIG,  29 
RESONANCE   LINE     SHAPE 


Fig.  27  is  a  composite  of  the  individual  line  shapes  shown  in 
Fig.  22.    As  previously  mentioned,  it  is  not  possible  to  separate  the 
absorption  and  dispersion  lines;  however,  since  a  pure  dispersion  line 
has  a  null  at  60=6jo    ,  the  peak  of  the  absorption  is  the  minimum  of  the 
trough  between  the  two  dispersion  peaks.    The  experimental  results  agree 
with  those  predicted  by  the  Bloch  susceptability  components.    This  leads 
to  the  conclusion  that  calculations  based  on  only  two  energy  levels  are 
sufficiently  accurate  to  represent  the  actual  system,  and  adds  a  measure 
of  confidence  to  the  extension  of  these  calculations  in  predicting  heading 
error . 
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The  graph  shown  in  Fig.  25  is  deceiving  in  that  it  suggests  that 
there  is  no  heading  error  dependence  on  the  phase  angle  between  the 
monitoring  light  and  coil  excitation.    For  sufficient  broadening,  the 
absorption  line  shape  is  practically  symmetrical  about  the  line  peak; 

o 

consequently,  introduction  of  a  phase  of  other  than  -90  causes  the 
frequency  of  oscillation  to  shift  to  one  side  of  the  line  maximum,    if 
this  shift, is,  say  +25  cycles/sec.  ,  then  by  changing  the  polarity  of 
circular  polarization  a  frequency  difference,  (Af-ss^  +  U-f+zs)  =  2Af 
is  obtained      Thus,  although  the  plotted  data,    2Af       ,  is  independent 
of  phase,  no  conclusions  can  be  extrapolated  for  heading  error. 
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EFFECT   OF  PHASE    SHIFT 
ON  HEADING   ERROR 


In  an  operating  magnetometer  the  phase  angle  between  the  light 

e 
if  fixed  at  -90;  therefore  little  valuable  information  is  lost  in  the 

experimental  procedure      Of  major  importance  is  the  fact  that  phase 

effects  can  be  discounted  in  the  results  of  the  primary  experiment  - 

heading  error  versus  orientation. 

It  has  been  assumed  that  the  principal  line  broadening  factor  is 

the  magnitude  of    H;     .    If  a  relatively  large  value  of    Hr      is  selected, 
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and  this  value  held  constant  as  the  magnetometer  is  rotated,  since 
H !  ■    HfSin-e-  ,  the  heading  error  should  vary  as  A-F=  Ap^^jHAsmdl I 

Fig.   30  is  a  plot  of  Af     versus  •©■     for    Hr  =  £"mv      .     For  large  values  of  •& 
it  was  necessary  to  extrapolate  the  value  of   Af    by  extending  a 
A^ipf  plot  to  an  r-f  value  of  5  mv. 

The  effect  of  the  optical  pumping  factor     \Pi+P*\  ,  on  heading 

error  is  shown  in  Fig.  31.    This  graph  is  the  result  of  varying  the 
orientation  while  holding   W{  constant.    Again,  the  incremental  values  of 
were  picked  off  of  the  Ia-c        Pl°t- 

Figs.  28  and  29  are  similar  to  the  curves  predicted  in  section  III, 
although  it  is  obvious  that  the  simple  relationships  stated  do  not  ade- 
quately represent  the  system. 

Another  factor  that  cannot  be  explained  at  this  time  is  the  relation 

between    H,     and  the  oscillating  frequency.    Referring  to  Fig.  26,  in  most 

cases  the  frequency  recorded  was  the  higher  of  the  frequencies  for  the 

two  polarizer  positions.    It  is  logical  that  as  the  r-f  drive  is  increased. 

the  line  is  broadened,  and  the  oscillating  frequency  should  decrease 

toward  f    .    If  the  theory  of  section  III  is  correct,  this  decrease  in  fre- 
o 

quency  should  be  one  half  the  decrease  in  ££p      .    This  is  not  the  case. 
For  -&=-50       ,  the  change  in  heading  error  is  10  cycles  while  the  change 
in  oscillating  frequency  is  7  2  cycles.     For  ■&■■=- 20°    ,  the  change  in  head- 
ing error  is  13  cycles  for  an  oscillating  frequency  change  of  11  cycles. 
This  would  seem  to  suggest  that,  since  for  larger  values  of  -&■  ,    H.  is 
more  nearly  equal     HT    ,  there  is  some  interaction  between  the  field  being 
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FIG-  30 

HEADING    ERROR    VER$u<S     ORIENTATION 
WITH     Hj    HELD  CONSTANT 
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measured  and  the  oscillating  r-f  field. 

Data  relating  to  the  dispersion  oscillation  does  not  plot  in  any 
significant  pattern.    Unfortunately,  the  primary  emphasis  of  these  experi- 
ments was  the  absorption  mode,  consequently,   the  data  taken  using  the 
dispersion  oscillation  was  sketchy.    The  general  trend,  however,  both 
in  the  data  recorded  and  in  many  observations  which  were  not  recorded , 
is  that  heading  error  is  considerably  less  for  the  dispersion  oscillation. 
This  goes  along  with  the  theory  that  heading  error  is  reduced  by  an 
increase  in  the  r-f  field. 
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VI.     CONCLUSIONS 

The  theoretical  development  and  experimental  results  both  show 
significant  effects  of  orientation  on  the  oscillating  frequency  of  the 
magnetometer,  and  although  only  the  system  using  Rb  85     was  investi- 
gated, the  theory  indicates  that  this  will  be  the  case  for  any  system 
using  an  absorption  media  having  nuclear  spin  other  than  zero. 

To  obtain  an  orientation  insensitive  instrument,  the  first  thought 
is  to  use  metastable  helium  which,  like  the  alkali  vapors,  is  charac- 
terized by  a  single  electron  in  the  outer  shell,  but  unlike  the  alkali 
vapors,  does  not  have  nuclear  spin.    This  solution  has  several  dis- 
advantages.   First,  the  high  potential  necessary  to  excite  the  helium 
to  the  metastable  state  gives  rise  to  additional  sources  of  error  such 
as  the  Stark  effect      Secondly,  for  a  magnetic  field  of  .5  gauss,  the 
Larmor  precession  frequency  of  metastable  helium  if  1.4  megacycles 
which  is  considerably  higher  than  the  frequency  response  of  present 
photo-detectors.    This  means  that  complex  detection  schemes  must  be 
devised  and  the  result  is  a  more  complicated  system.    There  are  other 
problems  connected  with  the  helium  magnetometer,  but  the  purpose 
here  is  just  to  show  that  a  switch  to  helium  to  eliminate  the  orientation 
effects  inherent  in  rubidium  is  probably  not  the  optimum  solution. 

Since  heading  error  in  the  rubidium  magnetometer  does  follow  a 
definite  pattern  which ,  with  a  little  more  experimentation ,  can  be  defined 
quite  precisely,  there  should  be  feasible  methods  by  which  it  can  be 
reduced  to  acceptable  limits . 
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One  of  these  methods  presently  being  studied  at  Varian  is  the 
use  of  multiple  cell  systems.    This  not  only  reduces  the  heading  error, 
but  also  eliminates  the  polar  and  equatorial  dead  zones  associated  with 
the  single  cell  unit. 

Another  possible  solution  is  to  build  an  orientation  sensitive  phase 
shifter  which  automatically  adjusts  the  oscillating  point  on  the  absorption 
line  to  compensate  for  any  frequency  shift  due  to  magnetometer  orientation. 
This  could  be  accomplished  by  a  sophisticated  switching  and  discriminator 
arrangement  that  would  sample  each  of  the  dispersion  frequencies  and 
produce  a  d-c  voltage  proportional  to  the  difference  frequency,  which, 
in  turn,  would  be  compared  to  a  d-c  voltage  proportional  to  the  r-f  input 
to  the  feedback  coil.    The  error  voltage  thus  produced  would  be  a  function 
of  the  magnetometer  orientation  since  the  separation  of  the  dispersion 
peaks  is  directly  related  to  the  component  of  HT  perpendicular  to     H Q     . 
The  error  voltage,  then,  could  be  applied  to  a  varicap  controlled  phase 
shifter  which  would  change  the  oscillator  frequency  to  the  correct  value. 

Still  another  solution  might  be  to  artificially  broaden  the  absorption 
line  by  adding  additional  impurities  to  the  rubidium  vapor.    This,  how- 
ever, may  not  be  a  wise  choice  in  that  it  is  equivalent  to  adding  noise 
to  the  system  and  would  deteriorate  the  performance  of  the  instrument. 

Further  experimentation  will  quite  possibly  produce  data  that  will 
suggest  better  ways  of  reducing  heading  error  than  those  just  discussed. 
With  this  in  mind,  it  may  be  beneficial  to  mention  some  additional  experi- 
ments that  the  writer  feels  are  essential  to  a  thorough  understanding  of  the 

problem . 
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1)  The  assumption  that  heading  error  is  symmetrical  about  the 
equatorial  dead  zone  must  be  verified.    This  may  be  done  by  mounting 
the  magnetometer  is  a  non-magnetic  cradle  and  the  oscillation  frequency 
observed  as  the  entire  system  is  rotated  360    .    This  is  not  a  simple 
experiment  in  that  the  most  carefully  designed  cradle  will  cause  some 
distortion  of  the  magnetic  field,  and  even  these  small  distortions  can 

-7 

be  significant  when  compared  to  the    \0        gauss  sensitivity  of  the 
instrument. 

2)  The  effect  of  phase  shift  on  the  oscillation  frequency  should 

be  determined.    Again,  this  is  not  a  simple  problem.    First  it  must  be 

o 
shown  that  the  assumptions  of  90  phase  shift  in  the  photo  cell  and 

negligible  phase  shift  in  the  r-f  coil  are  valid,  and  then  a  phase 

measuring  device  designed  that  will  measure  th^  phase  difference  between 

the  photo  detector  output,  (approx.    10      volts) ,  and  the  input  to  the  r-f 

coil,  (approx.  5    10       volts). 

3)  When  the  comparatively  large  effects  of  r-f  drive  and  phase 
shift  are  known  and  can  be  isolated,  then  smaller  effects,  such  as  those 
caused  by  light  intensity  and  absorption  cell  temperature,  can  be  investi- 
gated . 

This  final  section  was  written  under  the  assumption  that  it  is 
desirable  to  make  the  magnetometer  completely  orientation  insensitive. 
This  is  not  intended  to  convey  the  impression  that  the  rubidium  magne- 
tometer is  not  a  very  useful  instrument  regardless  of  heading  error. 
Actually  there  are  applications  where  useful  information,  such  as  the 
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rate  of  rotation  of  a  satellite,  can  be  obtained  by  making  use  of  the 
orientation  sensitive  characteristics  inherent  in  a  rubidium  device. 
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